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Abstract 


We present the discovery of a new Jovian-sized planet, TOI-3757 b, the lowest-density transiting planet known to orbit an 
M dwarf (MOV). This planet was discovered around a solar-metallicity M dwarf, using Transiting Exoplanet Survey 
Satellite photometry and confirmed with precise radial velocities from the Habitable-zone Planet Finder (HPF) and NEID. 
With a planetary radius of 12.0°$:3 R., and mass of 85.3188 M,, not only does this object add to the small sample of gas 
giants (~10) around M dwarfs, but also its low density (9 = 027 g cm *) provides an opportunity to test theories of 
planet formation. We present two hypotheses to explain its low density; first, we posit that the low metallicity of its stellar 
host (~0.3 dex lower than the median metallicity of M dwarfs hosting gas giants) could have played a role in the delayed 
formation of a solid core massive enough to initiate runaway accretion. Second, using the eccentricity estimate of 
0.14 + 0.06, we determine it is also plausible for tidal heating to at least partially be responsible for inflating the radius of 
TOI-3757b b. The low density and large scale height of TOI-3757 b makes it an excellent target for transmission 
spectroscopy studies of atmospheric escape and composition (transmission spectroscopy measurement of ~ 190). We use 
HPF to perform transmission spectroscopy of TOI-3757 b using the helium 10830 A line. Doing this, we place an upper 
limit of 6.9% (with 90% confidence) on the maximum depth of the absorption from the metastable transition of He at 
~10830 A, which can help constraint the atmospheric mass-loss rate in this energy-limited regime. 


Unified Astronomy Thesaurus concepts: Exoplanet astronomy (486); Hot Jupiters (753); Exoplanets (498); Radial 
velocity (1332); Exoplanet detection methods (489); Transits (1711) 
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1. Introduction 


Giant planets (R, > 4 Rg) should be intrinsically rare around 
M dwarfs according to planet formation models based on the 
core-accretion framework (Laughlin et al. 2004; Ida & Lin 
2005). In fact, recent simulations by Burn et al. (2021) find that 
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gas giants do not form for host stars <0.5M.. As the 
protoplanetary disk mass (and hence the amount of rocky 
material available) is correlated with the host-star mass 
(Andrews et al. 2013), disks around lower-mass stars such as 
M dwarfs should have less material available to form planetary 
cores. Laughlin et al. (2004) show that the low solid surface 
density of disks, coupled with the longer orbital timescales (due 
to lower stellar mass) makes it difficult to form these massive 
cores. If the cores become massive enough (>10 Mq; Pollack 
et al. 1996; Ida & Lin 2004) before the protoplanetary disk 
depletes, the process of runaway gas accretion is initiated, 
which is responsible for the formation of Jovian-type planets. 
In addition to a potential lack of disk mass around lower-mass 
stars, the difficulty of forming a Jovian-type planet around an 
M dwarf is further exacerbated by the shorter disk lifetimes in 
the UV-rich environment of M dwarfs, where the disks can 
rapidly deplete due to evaporation (Adams et al. 2004). 

Results from radial velocity (RV) searches (End et al. 2006; 
Johnson et al. 2007; Sabotta et al. 2021), transiting (Kovacs 
et al. 2013; Morton & Swift 2014; Obermeier et al. 2016), 
microlensing, and direct imaging (Gould et al. 2006; Montet 
et al. 2014) studies overall support this theory by constraining 
the occurrence rate of short-period gas giants around M dwarfs 
to ~1%-2%. New and ongoing giant planet confirmations from 
NASA’s Transiting Exoplanet Survey Satellite (TESS) mission 
(Ricker et al. 2014) will help refine these estimates further, with 
10 planets already added to this sample with TOI-1728b 
(Kanodia et al. 2020), TOI-1899b (Cafias et al. 2020), TOI- 
442b (Dreizler et al. 2020), TOI-674b (Murgas et al. 2021), 
TOI-532b (Kanodia et al. 2021), HATS-74Ab, HATS-75b 
(Jordan et al. 2021), and recently TOI-3629b and TOI-3714b 
(Cafias et al. 2022a). While the sample of M dwarf gas giants is 
small, several general trends appear to be emerging. One, the 
majority of M dwarf gas giants discovered by TESS are found 
orbiting early M-type stars, which should possess larger disks 
than their mid/late-M counterparts. Second, similar to the FGK 
orbiting hot Jupiters (Gonzalez 1997; Santos et al. 2001; 
Fischer & Valenti 2005; Ghezzi et al. 2010; Sousa et al. 2011), 
there is an apparent correlation between the occurrence of gas 
giants and host-star metallicity (Johnson & Apps 2009; 
Maldonado et al. 2019). However, pursuing this trend further 
with M dwarfs is currently hampered by their intrinsic faintness 
(in the optical) as well as the complexities with M dwarf 
metallicity determination (Passegger et al. 2022). 

An important distinction between M dwarf orbiting hot 
Jupiters and their FGK orbiting cousins is the large equilibrium 
temperature difference. As M dwarfs are significantly cooler, 
their hot Jupiters possess equilibrium temperatures <1000 K. It 
is therefore unlikely that M dwarf hot Jupiters experience the 
same inflation mechanism that “puffs” up the radius of the 
hotter hot Jupiters around FGK stars (e.g., ohmic dissipation; 
Batygin et al. 2011). That said, given that most M dwarf hot 
Jupiters fall within 0.05 au, they may experience some form of 
tidal heating assuming they are able to maintain a slightly 
eccentric orbit. Millholland et al. (2020) demonstrate that even 
a slightly eccentric orbit could explain the inflated radii of 
many low-density superpuffs, suggesting that these planets do 
not possess unusually large H/He atmospheres but hotter than 
expected interiors due to tidal forces. While this study has yet 
to be extended to M dwarf gas giants, due to the combination of 
lower stellar mass and radius (and hence a higher ratio of 
semimajor axis to stellar radius; a/R,.), these planets should 
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experience longer circularization timescales than their FGK 
counterparts. Therefore, it is possible that a subset of M dwarf 
hot Jupiters may currently be experiencing tidal heating in their 
interiors. 

In this manuscript we present the discovery of a low-density 
Jovian-sized planet with an inflated radius orbiting a solar- 
metallicity M dwarf—TOI-3757 in the constellation of Auriga. 
We use a combination of photometry from TESS and ground- 
based instruments (RBO), high-contrast speckle imaging 
(NESSD, and precision RVs from the Habitable-zone Planet 
Finder (HPF; Mahadevan et al. 2012, 2014) and NEID 
(Halverson et al. 2016; Schwab et al. 2016) spectrographs. 
We also use HPF to observe the planet during its transit and 
perform transmission spectroscopy to place upper limits on 
absorption in He 10830A. In Section 2 we detail these 
observations which are used to characterize the system. In 
Section 3 we describe the methodology followed to derive the 
stellar parameters, while Section 4.1 details the data analysis, 
including the joint fitting of the photometry and RVs, and also 
the upper limits we place on He _ 10830 A absorption 
(Section 4.2). In Section 5 we place TOI-3757 b in context 
of other planets around M dwarfs and also hypothesize 
different mechanisms that can be responsible for its low 
density, while finally summarizing our findings in Section 6. 


2. Observations 
2.1. TESS 


TOI-3757 (TIC-445751830, Two Micron All Sky Survey 
(2QMASS) J06040089+5501126, Gaia Early Data Release 
(EDR) 3 996878131494639488, UCAC4 726-038940) is an 
MO dwarf observed by TESS in Sector 19 in Camera 2 from 
2019 November 27 to 2019 December 24 at ~30 minute 
cadence (Figure 1). The planet candidate was identified using 
the Quick Look Pipeline (QLP) algorithm developed by Huang 
et al. (2020), under the “faint-star search” (Kunimoto et al. 
2021) with a period of ~3.43 days. 

We extract the photometry from the TESS full-frame images 
using eleanor (Feinstein et al. 2019), which uses the 
TESScut* service to obtain a cutout of 31 x 31 pixels from 
the calibrated full-frame images centered on TOI-3757. The 
TESS light curve is derived from the CORR_FLUX values, in 
which eleanor uses linear regression with pixel position, 
measured background, and time to remove signals correlated 
with these parameters. The final aperture is a 2 x | rectangle 
centered on TOI-3757 and was selected by minimizing the 
combined differential photometric precision (CDPP) after the 
data were binned in | hour timescales. The CDPP is formally 
the rms of the photometric noise on transit timescales and was 
originally defined for Kepler (Jenkins et al. 2010). We obtain a 
CDPP of 2730 ppm for the TESS photometry. 


2.2. Ground -based Photometric Follow Up 
2.2.1. Red Buttes Observatory 


We observed a transit of TOI-3757 b on the night of 2021 
November 17 using the 0.6m telescope at the Red Buttes 
Observatory (RBO) in Wyoming (Kasper et al. 2016). The 
telescope is a f/8.43 Ritchey—Chrétien Cassegrain constructed 


3° https: //mast.stsci.edu/tesscut/ 
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Figure 1. Long-cadence (30 minute) time-series TESS eleanor photometry from Sector 19, along with a stellar rotation GP kernel (RotationTerm from 
celerite2) in green. The detrended photometry is shown in the bottom panel, with the TOI-3757 b transits overlaid in blue. 


by DFM Engineering, Inc, and it is currently equipped with an 
Apogee Alta F16 camera. 

The target rose from an airmass of 1.20 at the start of the 
observations to a minimum airmass of 1.02 and then set to an 
airmass of 1.04 at the end of the observations. We performed 
defocussed observations using the Bessell I filter (Bessell 1990) 
exposure times of 240s. In the 2 x 2 binning mode, the 0.6m 
at RBO has a gain of 1.39 e/ADU, a plate scale of 073, and a 
readout time of approximately 2.4 s. 

For the final reduction, we selected a photometric aperture of 
13 pixels (95) with a sky annulus of inner and outer radius of 
23 pixels (16"8) and 35 pixels (255), respectively. We obtain 
an rms precision of ~1750 ppm on the RBO photometry, after 
subtracting the transit model (Figure 6). 


2.3. NN-Explore Exoplanet Stellar Speckle Imager at WIYN 


We observed TOI-3757 on the night of 2021 December 21 
using the NN-Explore Exoplanet Stellar Speckle Imager 
(NESSI) on the WIYN 3.5m telescope at Kitt Peak National 
Observatory to search for faint background stars and nearby 
stellar companions. A 9 minute sequence of 40 ms diffraction- 
limited exposures was collected using the Sloan r’ filter on 
NESSI. The speckle images were reconstructed following the 
procedures described in Howell et al. (2011). No stellar sources 
were detected down to a magnitude limit of Ar’=4.0 at 
separations >0”2, as shown in Figure 2. 


2.4. Radial Velocity Follow Up 
2.4.1. HPF 


TOI-3757 was observed using HPF (Mahadevan et al. 
2012, 2014) starting 2021 September 1. HPF is a near-infrared 
(8080-12,780 A), high-resolution fiber-fed (Kanodia et al. 
2018) precision RV spectrograph with exceptional environ- 
mental stability (Stefansson et al. 2016) located at the 10 m 
Hobby-Eberly Telescope (HET) at McDonald Observatory, 
Texas. HET is a fixed-altitude telescope with a roving pupil 
design and is fully queue-scheduled, where all the observations are 
executed by the HET resident astronomers (Shetrone et al. 2007). 
Using the algorithms described in the package HxRGproc (Ninan 
et al. 2018), we correct for bias, nonlinearity, cosmic rays, and 
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Figure 2. 5c contrast curve for TOI-3757 observed from NESSI in the Sloan r’ 
filter showing no bright companions within 1/’2 from the host star. The r’ 
image is shown as an inset 1” across. 


calculate the slope/flux and variance images from the raw HPF 
data. While HPF has the capability for simultaneous calibration 
using a NIR Laser Frequency Comb (LFC; Metcalf et al. 2019), 
due to the faintness of our target we chose to avoid simultaneous 
calibration to minimize the impact of scattered calibrator light in 
the science target spectra. Instead, we obtain a wavelength solution 
for the target exposures by interpolating the wavelength solution 
from other LFC exposures on the night of the observations. This 
helps correct for the well-calibrated instrument drift (Stefansson 
et al. 2020). This method has been shown to enable precise 
wavelength calibration and drift correction with a precision of ~30 
cms! per observation, a value much smaller than our estimated 
per observation RV uncertainty (instrumental + photon noise) for 
this object of 34 ms! (in 969 s exposures and 23 ms! in binned 
30 minute exposures). 

We follow the method described in Stefansson et al. (2020) 
to derive the RVs, by using a modified version of the 
SpEctrum Radial Velocity AnaLyser pipeline (SER- 
VAL; Zechmeister et al. 2018). SERVAL uses the template- 
matching method (e.g., Anglada-Escudé & Butler 2012), where 
it creates a master template from the target star observations, 
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Table 1 
RVs (Binned in ~30 minute Exposures) of TOI-3757 

BJDrpp RV oO Instrument 
(d) ms! ms! 

2459458.98300 —28.29 22.93 HPF 
2459466.97200 42.47 24.02 HPF 
2459472.95600 —11.67 23.40 HPF 
2459489.92300 —62.93 21.99 HPF 
2459503.87500 —36.31 28.99 HPF 
2459506.86400 —71.04 23.27 HPF 
2459507.87000 —22.90 20.91 HPF 
2459509.85300 —50.03 23.73 HPF 
24595 12.83800 —75.03 22.47 HPF 
24595 13.83100 —109.27 28.70 HPF 
24595 16.84000 —96.78 18.73 HPF 
24595 17.83600 —52.47 19.65 HPF 
2459532.80200 —2.53 26.69 HPF 
2459533.78400 —67.68 28.71 HPF 
245955 1.75800 —139.99 25.00 HPF 
2459571.69600 —102.35 25.02 HPF 
2459505.96900 —0.50 17.21 NEID 
2459523.93279 —36.64 6.78 NEID 
2459528.80595 66.52 10.79 NEID 
2459531.81514 27.19 11.85 NEID 
2459532.81927 46.97 10.55 NEID 
2459534.03339 —43.36 11.05 NEID 
2459546.785 16 55.86 12.63 NEID 
2459554.99677 —21.28 Tat NEID 
2459582.88049 19.91 12.05 NEID 
2459585.89929 —20.71 8.26 NEID 
2459589.89730 —1.91 11.30 NEID 


Note. We include this table in a machine-readable format along with the 
manuscript. 


(This table is available in its entirety in machine-readable form.) 


and determines the Doppler shift for each individual observa- 
tion by minimizing the x? statistic. The master template is 
created using all of the HPF observations for TOI-3757, after 
masking out the telluric and sky emission lines. The telluric 
regions are identified by a synthetic telluric-line mask 
generated from telfit (Gullikson et al. 2014), a Python 
wrapper to the line-by-line radiative transfer model package 
(Clough et al. 2005). We use barycorrpy to perform the 
barycentric correction on the individual spectra, which is the 
Python implementation (Kanodia & Wright 2018) of the 
algorithms from Wright & Eastman (2014). 

We obtained a total of 25 visits on this target between 2021 
September | and 2021 December 24, of which 9 had to be 
excluded from further analysis due to poor weather conditions 
during the observations (seeing and sky transparency).°*! Each 
visit was divided into two exposures of 969s each, where the 
median S/N of each HPF exposure was 47 per pixel at 
1070 nm. The individual exposures were then binned after 
weighting based on SNR, with the final binned RVs being 
listed in Table 1. 

Additionally, we observed TOI-3757 with HPF on 2021 
December 12, during the transit of planet b. This observation 
consisted of 10x individual exposures of 649s each and was 


3! The rejected visits had an unbinned median S/N and RV uncertainty of ~32 
and ~53 msl, compared to 47 and 34 ms |, respectively, for the visits 
included in the analysis. 
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used to place constraints on atmospheric escape using the 
infrared atomic transitions of helium as a tracer (discussed 
further in Section 4.2). Out of an abundance of caution, we do 
not include this visit in our RV analysis, to avoid potential 
systematics from the Rossiter-McLaughlin effect (Rossi- 
ter 1924; McLaughlin 1924; Triaud 2018). 


2.4.2. NEID 


We also observed TOI-3757 using NEID, a new ultra-precise 
(Halverson et al. 2016), environmentally stabilized (Robertson 
et al. 2019) spectrograph at the WIYN 3.5 m telescope at Kitt 
Peak National Observatory. NEID is a_ high-resolution 
spectrograph (R ~ 110, 000) with an extended red wavelength 
coverage (380-930 nm; Schwab et al. 2016); it has a fiber-feed 
system similar to HPF (Kanodia et al. 2018) with three fibers— 
science, sky, and simultaneous calibration. For these observa- 
tions we use the NEID high-resolution (HR) mode** with 
resolution R ~ 110, 000. Between 2021 November 1 and 2022 
January 10, we obtained 12 visits on TOI-3757 with NEID, of 
which we exclude the visit from 2021 November 2 with a S/N 
of 2.3 at 850 nm due to patchy clouds during the observation. 
Of the visits included for analysis, each consisted of an 1800s 
exposure with a median S/N per ID extracted pixel of 9.3 at 
850 nm. 

The NEID data were reduced using the NEID data reduction 
pipeline** (DRP), and the Level-2 1D extracted spectra were 
retrieved from the NEID archive.** We used a modified version 
of the SERVAL template-matching algorithm to obtain the 
NEID RVs (Stefansson et al. 2022), similar to that for HPF. 
The NEID RVs presented here were calculated using orders 
spanning 4560-8960 A (order indices 40 to 104), across the 
central 7000 pixels which masks out the low S/N spectra 
outside the free spectral range. We also mask out the telluric 
and sky emission lines, which were identified following an 
identical procedure to our HPF processing (Section 2.4.1), and 
obtain the barycentric velocities using Kanodia & Wright 
(2018). The final NEID RVs are listed in Table 1. 

The early spectral type (MO) of TOI-3757 leads to higher RV 
information in the optical than in the NIR (Bouchy et al. 2001). 
Coupled with negligible rotational broadening (v sini << 2kms_'), 
this leads to NEID RVs (R~ 110, 000) that are much more 
precise than HPF (R~ 55, 000; Table 1, Figure 3). 


3. Stellar Parameters 
3.1. Spectroscopic Parameters with HPF-SpecMatch 


We use the HPF-SpecMatch* package (Stefansson et al. 
2020) to empirically determine stellar parameters from HPF 
spectra using the template-matching method based on Yee et al. 
(2017). 

The HPF-SpecMatch algorithm employs a two-step process, 
which uses the a metric twice to find the library stars that best fit 
the target spectrum. The HPF spectral library contains 166 stars 
and spans the following parameter space: 2700 K < T, < 6000 K, 
4.3 < logg, < 5.3, and —0.5 < [Fe/H] < 0.5. 

In the first step, each stellar library spectrum is compared 
with the target spectrum using the y” metric and ranked from 


3? Instead of the high-efficiency (HE) mode with resolution R ~ 70, 000. 
33 https: / /neid.ipac.caltech.edu /docs/NEID-DRP/ 

34 https: / /neid.ipac.caltech.edu/ 

= https: //gummiks.github.io/hpfspecmatch / 
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Figure 3. Time series of RV observations of TOI-3757 with HPF (green) and NEID (red). The best-fitting model derived from the joint fit to the photometry and RVs 
is plotted in blue, including the 16%—84% confidence interval in lighter blue. The bottom panel shows the residuals after subtracting the model. 


best to worst-fitting library star to the target spectrum. In the 
second step, only the top five best-fit ranked library stars are 
used. The y~ metric is applied to assign scaling constants to 
each of the five best-fit library stars and create a composite 
spectrum that fits the target spectrum even more closely. These 
scaling constants are also used to determine a weighted average 
for the precise parameter estimates of effective temperature 
Tere, surface gravity (log g), and metallicity ({[Fe/H]) of the 
target star. 

The errors we adopt for the spectroscopic parameters are 
determined via a leave-one-out cross validation, broadly 
following the methodology in Stefansson et al. (2020). In this 
process, a library star of interest is removed from the rest of the 
stellar library pool. Then, the HPF-SpecMatch algorithm is 
run to estimate its stellar parameters independently of its true 
parameter. The difference between these calculated parameters 
and the true parameter values are noted. The spectral matching is 
performed on HPF order index 5 (8534-8645 A) for TOI-3757 
because this order has negligible telluric contamination. The 
resolution limit of HPF places a constraint of v sini < 2kms7! 
for TOI-3757. Table 2 presents the derived spectroscopic 
parameters with their uncertainties. 


3.2. Model-dependent Stellar Parameters 


We derive model-dependent stellar parameters by modeling 
the spectral energy distribution (SED) for TOI-3757 using the 
EXOFASTv2 analysis package (Eastman et al. 2013). The SED 
fit uses the precomputed bolometric corrections*® in log g,, Tete, 
[Fe/H], and Ay from the Modules for Experiments in Stellar 
Astrophysics Isochrones and Stellar Tracks (MIST) model 
grids (Dotter 2016; Choi et al. 2016). 

We place Gaussian priors on the (i) broadband photometry 
listed in Table 2, (ii) the spectroscopic stellar parameters 
derived with HPF-SpecMatch, and (iii) the geometric 
distance calculated from Bailer-Jones et al. (2021). We set 


36 http: //waps.cfa.harvard.edu/MIST /model_grids.html#bolometric 


the upper limit of the visual extinction to the estimate from 
Green et al. (2019) calculated at the distance determined by 
Bailer-Jones et al. (2021). The extinction from Green et al. 
(2019) is converted to a visual magnitude extinction using the 
R,=3.1 reddening law from Fitzpatrick (1999). Table 2 
contains the stellar priors and derived stellar parameters with 
their uncertainties. 


3.3. Galactic Kinematics 


We use the systemic velocity derived from HPF and proper 
motion from GAIA EDR3 to calculate the UVW velocities in 
the barycentric frame using GALPY (Bovy 2015).°” We also 
provide the UVW velocities in the local standard of rest using 
the offsets from Sch6nrich et al. (2010). Using the BANYAN 
tool (Gagné et al. 2018), we classify TOI-3757 as a field star in 
the thin disk with very high probability (>99%; Bensby et al. 

2014). 


3.4. Rotation Period Estimates 


We run a generalized Lomb-Scargle (GLS) periodogram 
(Lomb 1976; Scargle 1982; Zechmeister & Kiirster 2009) on 
the TESS photometry (after masking the transits on TOI-3757 
b) using its astropy implementation and find a significant 
peak (<0.1% False Alarm Probability) at ~7 days (Figure 4). 
This is consistent with the results from our Gaussian process 
(GP) stellar rotation kernel** applied to the TESS photometry, 
which suggest a stellar rotation period of 6.9')3 days. This 
kernel consists of two simple harmonic oscillator terms—one at 
the rotation period, with the second one at half the period. 

On further inspection using eleanor we find a similar 
significant peak at ~7 days for the photometry in the majority 
of the adjoining pixels in a 6 x 6 grid centered on the centroid 


37 With U toward the Galactic center, V toward the direction of Galactic spin, 
and W toward the North Galactic Pole (Johnson & Soderblom 1987). 

38 RotationTerm implemented in celerite2 (Foreman-Mackey et al. 
2017; Foreman-Mackey 2018). 
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Table 2 

Summary of Stellar Parameters for TOI-3757 
Parameter Description Value Reference 
Main identifiers: 
TOI TESS Object of Interest 3757 TESS mission 
TIC TESS Input Catalogue 445751830 Stassun 
2MASS a J06040089+-5501126 2MASS 
Gaia EDR3 vee 996878 131494639488 Gaia EDR3 
Equatorial Coordinates, Proper Motion, and Spectral Type: 
232016 R.A. (RA) 06:04:00.87 Gaia EDR3 
632016 decl. (Dec) 55:01:11.90 Gaia EDR3 
La Proper motion (RA, mas/yr) —9.03 + 0.02 Gaia EDR3 
Lbs Proper motion (Dec, mas/yr) —43.13 + 0.02 Gaia EDR3 
d Distance in pc 177.4 + 0.7 Bailer-Jones 
Ay max Maximum visual extinction 0.26 Green 
Optical and Near-infrared Magnitudes: 
B Johnson B mag 16.2 +0.2 APASS 
V Johnson V mag 14.8+0.1 APASS 
g! Sloan g’ mag 15.5+0.1 APASS 
r! Sloan r’ mag 14.2+0.1 APASS 
i! Sloan i’ mag 13.5+0.1 APASS 
J J mag 12.00 + 0.03 2MASS 
H H mag 11.31 + 0.03 2MASS 
K, K, mag 11.15 + 0.02 2MASS 
wil WISE! mag 11.06 + 0.02 WISE 
w2 WISE2 mag 11.10 + 0.02 WISE 
Ww3 WISE3 mag 11.0+0.1 WISE 
Spectroscopic Parameters": 
Tete Effective temperature in K 3913 + 56 This work 
[Fe/H] Metallicity in dex 0.0 + 0.20 This work 
log(g) Surface gravity in cgs units 4.68 + 0.04 This work 
Model-dependent Stellar SED and Isochrone fit Parameters”: 
My. Mass in M, 0.64 + 0.02 This work 
Ry. Radius in Re 0.62 + 0.01 This work 
Ly Luminosity in L 0.087 + 0.003 This work 
Px Density in g cm 3.7+0.2 This work 
Age Age in Gyrs 71445 This work 
A, Visual extinction in mag 0.067+0-078 This work 
Other Stellar Parameters: 
v Sinix Rotational velocity in km s~! <2kms! This work 
ARV “Absolute” radial velocity in km s~! 21.86 + 0.04 This work 
U, V, W Galactic velocities in km s~! 36.12+0.07, —17.8140.10, —15.82+0.09 This work 
U, V, W Galactic velocities (LSR) in km s 25.02+0.85, —5.5740.69, —8.56+0.61 This work 


Note. References are: Stassun (Stassun et al. 2018), 2MASS (Cutri et al. 2003), Gaia EDR3 (Collaboration et al. 2021), Bailer-Jones (Bailer-Jones et al. 2018), Green 
(Green et al. 2019), American Association of Variable Star Observers Photometric All Sky Survey (APASS; Henden et al. 2018), Wide-field Infrared Survey Explorer 


(WISE; Wright et al. 2010). 


* Derived using the HPF spectral matching algorithm from Stefansson et al. (2020). 


> EXOFASTv2 derived values using MIST isochrones with the Gaia parallax and spectroscopic parameters in a) as priors. 
© The barycentric UVW velocities are converted into local standard of rest (LSR) velocities using the constants from Schénrich et al. (2010). 


for TOI-3757. This suggests that the periodic signal seen in the 
eleanor reduction of the TESS FFI photometry is not 
astrophysical in origin, and TOI-3757 does not have significant 
rotational modulation. This is further corroborated by the lack 
of a detectable rotational broadening signal in the HPF spectra, 
using which we can place a limit of vsini<2kms_' on the 
host star.* 

We also access publicly available data from the Zwicky 
Transient Facility (ZTF; Masci et al. 2019) and ASAS-SN 
(Kochanek et al. 2017) for this target to perform a GLS 
periodogram analysis, and do not detect any rotation signal 


° The corresponding equatorial velocity for a ~7 day rotation period would 
be ~4.5 kms’. 


present in the photometry. The photometry spans ~800 days for 
ZTF and ~1000 days for ASAS-SN. We also analyze the Ca 
infrared triplet from the HPF spectra, and Ha from the NEID 
spectra, but do not find any periodic signals in the time series. 


3.5. Ruling out Stellar Companions 
3.5.1. Unresolved Stellar Companions 


We try to place constraints on the presence of unresolved 
stellar companions using HPF spectra, Gaia astrometry and 
the RVs: 


1. Constraints from HPF spectra: We follow the procedure 
outlined in Kanodia et al. (2020), to place limits on any 
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Figure 4. GLS periodogram for the different data sets. The horizontal gray 
lines represent the 0.1%, 1%, and 10% FAP values. The red vertical line 
depicts the orbital period of planet b, whereas the gray line marks its 1 day 
alias. 


spatially unresolved stellar companion to TOI-3757 using 
the HPF spectra to quantify the lack of flux from a 
secondary object. We combine the spectra from a single 
epoch to obtain a higher S/N template for comparison 
and then model the test spectra (TOI-3757) as a linear 
combination of a primary M dwarf (BD+29_2279) and a 
secondary companions (GJ 251, GJ 1156 and VB-10). 
The flux ratio between the secondary and primary star, F, 
is calculated as: 


Sobs = AC = X) Sprimary oF (x) Ssecondary)s (1) 


F=— (2) 


l—x’ 


where Sops is the observed spectrum, Sprimary 18 the 
primary spectrum, Ssecondary represents the secondary 
spectrum, and A is the normalization constant. For a 
given primary and secondary template, we (i) perform a 
x minimization to shift the secondary spectrum in 
velocity space, (ii) add this shifted secondary spectrum to 
the primary, and (iii) fit for the value of x (and A) that best 
fits the observed spectrum. We perform this for a range of 
spectral types for the secondary from M3 to M8 spanning 
velocity offsets of +150 km s |. We place a conservative 
upper limit for a secondary companion of flux ratio <0.15 
or Amag ~ 2.1 for |Av| > 5 km s~', using HPF order index 
17 spanning 10450-10580 A. The lower limit coincides 
with HPF’s spectral resolution (R ~ 55, 000 = 5.5 kms_’). 
At lower velocity offsets, the degeneracy between the 
primary and secondary spectra prevents any meaningful flux 
ratio constraints. 

2. Constraints from Gaia astrometry: GAIA EDR3 (Lindegren 
et al. 2021) provides an additional astrometric constraint 
on the presence of unresolved bound companions using 
the renormalized unit weight error (RUWE) metric. 
RUWE is sensitive to the change in the position of the 
primary target due to reflex motion caused by unresolved 
bound companions. For the single-star astrometric solu- 
tion in use for GAIA EDR3, this astrometric motion of 
the primary star around the center of mass would manifest 
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as noise (Kervella et al. 2019), especially for orbital 
periods much shorter than the observing baseline for 
GAIA EDR3 (~34 months). The commonly accepted 
threshold in literature for this is RUWE= 1.4, which 
correlates with the presence of a bound stellar companion 
in recent studies of stellar binaries (Penoyre et al. 2020; 
Belokurov et al. 2020; Gandhi et al. 2021). For TOI-3757, 
GAIA EDR3 reports a RUWE of ~1.1, which is in 
agreement with a single-star astrometric solution. 

3. Constraints from RVs: We use the joint fit of the 
photometry and radial velocity to estimate the planetary 
and system properties (Section 4.1). We also include a 
linear RV trend in the orbital solution while fitting the 
RVs, and note this to be consistent with 0, where the RV 
trend w— 0,23 (ms! yr~'). The residuals to this fit 
(shown in Figure 3), are also analyzed with a GLS 
periodogram and show no significant signals (Figure 4), 
indicating the absence of any long-period bound 
companions over our observing baselines (~100 days). 


3.5.2. Resolved Stellar Companions 


Figure 5 presents a comparison of the region contained in 
the 11 x 11 pixel footprint from Sector 19 using a Palomar 
Observatory Sky Survey (POSS-1; Harrington 1952; Minkowski 
& Abell 1963) image from 1951 and a ZTF (Masci et al. 2019) 
image from 2019. The POSS-1 plate images were taken with 
Eastman 103a-E spectroscopic plates in conjunction with a No. 
160 red plexiglass filter with a bandpass between 6000 and 
6700 A, and have a limiting magnitude of R~ 19 (Harring- 
ton 1952). TOI-3757 has low proper motion, and the change in 
coordinates between the two epochs is negligible. There are no 
bright targets with A Grp < 3 present in the 2x 1 TESS 
aperture. There are a few targets with A Gpp < 4 that may dilute 
the TESS transit. We use the ground-based photometry to 
estimate the dilution term for the TESS photometry. Addition- 
ally, using the NESSI observations we are able to rule out stellar 
sources with a Ar’ < 4.0 at separations >0/2. 


4. Data Analysis 
4.1. Joint Fitting of Photometry and RVs 


We perform a joint fit of the photometry (TESS + ground- 
based sources) and the RVs (HPF + NEID) using the Python 
package exoplanet, which builds upon PyMC3, the 
Hamiltonian Monte Carlo (HMC) package (Salvatier et al. 
2016). The exoplanet package uses starry (Luger et al. 
2019; Agol et al. 2020) to model the planetary transits, using 
the analytical transit models from Mandel & Agol (2002), and a 
quadratic limb-darkening law. These limb-darkening priors are 
implemented in exoplanet using the reparameterization 
from Kipping (2013) for uninformative sampling. We fit each 
phased transit (Figure 6) with separate limb-darkening coefficients. 
Our likelihood function for the TESS photometry includes the GP 
kernel to model the stellar rotation signal. To account for the long- 
cadence photometry from TESS, exoplanet (Foreman-Mackey 
et al. 2021a) oversamples the time series while evaluating the 
model. 

We model the RVs using a standard Keplerian model, where 
we let the eccentricity float. We include a separate RV offset 
for each instrument (HPF and NEID), along with a common 
linear RV trend to account for long-term drifts. For the joint 
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Figure 5. Panel (a) overlays an 11 x11 pixel footprint from TESS Sector 19 (blue grid) on a POSS-I red image from 1953.1. TOI-3757 has a small proper motion as 
can be seen while comparing panels (a) and (b). The TESS aperture is outlined in red, and we highlight TOI-3757 with a star. No bright targets are present inside the 


TESS aperture with A Gpp < 3. Panel (b) is similar to panel (a) but with a background image from ZTF zr (5600-7316 A) from 2019 (Masci et al. 2019). 
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Figure 6. Photometric observations for TOI-3757 b. (a) the TESS-phased plot shows the detrended light curve phase-folded to the best-fit orbital period; (b) ground- 
based observations for TOI-3757 b from RBO. The raw photometry is shown in gray, with the best-fit transit solution in blue, along with the 1o confidence interval 


shown in lighter blue. The data behind Figure 6(b) are available. 


(The data used to create this figure are available.) 


modeling of the photometry + RVs, we also include a simple 
white-noise model in the form of a jitter term that is added in 
quadrature to the measurement errors from each data set. 
Using scipy.optimize, we find the initial maximum 
a posteriori (MAP) parameter estimates, which uses the default 
Broyden—Fletcher—Goldfarb-Shanno (BFGS) algorithm (Broy- 
den 1970; Fletcher 1970; Goldfarb 1970; Shanno 1970). These 
parameter estimates are then used as the initial conditions for 
parameter estimation using “No U-Turn Sampling” (NUTS; 
Hoffman & Gelman 2014), implemented for the HMC sampler 


PyMC3, where we check for convergence using the Gelman— 
Rubin statistic (R < 1.1; Ford 2006). 

The final derived planet parameters from the joint fit are 
included in Table 3, with the phased RVs shown in Figure 7. 
We obtain a ~~ 10a mass for TOI-3757 b of 85.3183 Mz, and a 
radius of 12.0193 Ro. 

From the joint RV + photometry fit we obtain an 
eccentricity of 0.14+0.06, which is consistent with the 
eccentricity estimate derived from the RBO photometry using 

the photoeccentric effect (Dawson & Johnson 2012). Here, we 
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Parameter 


Orbital Parameters: 
Orbital Period... 
Eccentricity... 


Argument of Periastron... 
Semi-amplitude Velocity... 


Systemic Velocity”... 


RV trend... 
RV jitter... 


Transit Parameters: 
Transit Midpoint ... 
Scaled Radius... 


Scaled Semimajor Axis... 


Orbital Inclination... 
Transit Duration... 
Photometric Jitter® ... 


Limb Darkening® ... 


Table 3 
Derived Parameters for the TOI-3757 System 


Units 


P (days) ... 
e@... 

w (degrees) ... 
K(ms_})... 
‘apr (ms *)... 
Yew (ms—')... 
dv/dt (ms! yr~') 
OHPE (m s 7)... 
Onew (ms-')... 


Tc (BJDzppz)..-- 
R, [Ry «+ 
a/Ry ... 

i (degrees)... 

T,4 (days)... 
Cress (ppm)... 
ORBo20211117 (ppm)... 
U\,TESS> 42,TESS -:- 


41,RBO20211117> 42,RBO20211117 +++ 
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Value* 


3.438753+0.000004 
0.14+0.06 
130423 
49.24+5.07 
—37.914789 
17.5844.15 
—0.207428 
17.71.5595" 
5.821366 


2458838.77 1484 0.00092 
0.1769 40.0056 
13.2618 
86.761033 
0.0800+0.0038 
2569+ 8 
2953+333 
0.397041, 0.102038 
0.594043, 0.071043 


Dilution’... Dress --- 0.994+0.053 
Planetary Parameters: 

Mass... M, (Mg)... 85.3788 
Radius... R, (Re). 12.0734 
Density... Pp (gem)... 0.274083 
Semimajor Axis... a (au) ... 0.03845+0.00043 
Average Incident Flux’... (F) (10° W/m”)... 0.75-+0.05 
Planetary Insolation S (Sg)... 55.4+£3.8 
Equilibrium Temperature® ... Teg (K)..- 759+ 13 


Note. 
*“ The reported values refer to the 16%-50%-84% percentiles of the posteriors. 


> In addition to the absolute RV from Table 2. 
© Jitter (per observation) added in quadrature to photometric instrument error. 
4 Where uy + uz < 1, and u, > 0 according to Kipping (2013). 


* Dilution due to the presence of background stars in the TESS aperture, not accounted for in the eleanor flux. 
We use a solar flux constant = 1360.8 W m~ to convert insolation to incident flux. 
© We assume the planet to be a blackbody with zero albedo and perfect energy redistribution to estimate the equilibrium temperature. 


compare the stellar density obtained from isochrones and SED 
fitting in Section 3, to that from the measured transit duration, 
which is affected by the eccentricity of this orbit. Using this 
comparison, we measure the eccentricity of the orbit from a 
joint fit of the TESS and RBO photometry to be 0. 127, as This 
is consistent (to within 1a) with the eccentricity derived from 
the joint RV + photometry fit (Table 3). 


4.2. Upper Limit on Helium 10830 A absorption 


The low bulk density of this planet makes it a promising 
candidate to detect mass loss via atmosphere escape. We 
observed TOI-3757 b during its transit with HPF (described in 
Section, 2.4.1), to constrain the absorption in the He 
10830 A triplet due to the planetary exosphere. Not only does 
HET’s fixed-altitude design drastically reduces the number of 
transits observable with HPF, but also the stellar barycentric 
velocity shifts the helium triplet feature against the bright 
hydroxyl sky emission lines. Considering these restrictions, we 
observed the best available transit of TOI-3757 b with HPF in 
2021 December. 


The individual spectra were sky subtracted using the 
simultaneous sky spectra obtained via an adjacent sky fiber 
(Kanodia et al. 2018). Conservatively, we inflate the error bars 
in the pixels that are corrected for sky emission lines to avoid 
potential systematics. Telluric correction was not performed on 
this spectra, as the wavelength region of interest for He 
10830 A detection was well outside the telluric absorption 
lines. We considered all spectra taken 1.9 hr (one full transit 
duration) away from transit midpoints in the 2021 November 
and December for the out-of-transit spectra. The individual out- 
of-transit spectra were then combined by weighted averaging to 
obtain a high signal-to-noise ratio template spectrum of the star. 
Individual in-transit spectra were then divided by this template to 
obtain a set of transmission (or ratio) spectra. These ratio spectra 
were then aligned to the planet’s rest frame (based on the orbit 
model) and combined by weighted average to obtain the final 
transmission spectrum shown in Figure 8. No statistically 
significant He 10830 A absorption signal was detected. 

For calculating an upper limit, we used a Gaussian 
absorption line model with an FWHM width of 0.89 A. This 
is the typical width of the reported helium detections in the 
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Figure 7. HPF + NEID RV observations phase-folded on the best-fit orbital 
period from the joint fit from Section 4.1. The best-fit model is shown in the 
solid line, whereas the lo confidence intervals are shown in lighter blue. 


6.9% upperlimit 
In-transit/Out-of-transit 


Normalized In-transit/Out-of-transit spectrum 


Wavelength [A] (in planet rest frame) 


Figure 8. Ratio of the in-transit spectra and out-of-transit TOI-3757 spectra. 
The blue curves are the three individual ratio spectra from the transit epoch, 
whereas the black curve is the weighted average of the three. The x-axis shows 
a vacuum wavelength in the planet’s rest frame mid-transit. The rest vacuum 
wavelengths of the He 10830 A, triplet lines in planet’s rest frame are marked 
by dashed vertical orange lines. We do not detect any significant absorption in 
the planetary spectra at these wavelengths. The results of our MCMC fit of the 
strongest doublet lines in the He 10830 A triplet using a Gaussian model of 
width 0.89 A are shown by the red curves in the lower panel, and the 6.9% 
upper limit is shown by the dashed red curve overlaid on the MCMC results. 


literature across different planets.*° For a circular orbit, we 
obtain an upper limit on the maximum depth of the line as 6.9% 
(with 90% confidence; see Figure 8). If we consider an 
eccentricity of 0.14 and argument of periastron to be 130 
degrees, the net radial velocity of the planet during the transit is 
redshifted by 11.2 kms “ in comparison to the circular orbit. 
This places the expected planetary signal on top of a bright 
hyroxyl sky emission line. Therefore, an eccentric orbit 
ephemeris precludes us from placing any meaningful upper 
limit. With this caveat, while our upper limit shows the signal 
is weaker than some of the strongest known signal in 
exoplanets (Allart et al. 2019; Kirk et al. 2020), we encourage 


40 EWHM of He 10830 Adetections in WASP 69b: 0.72 A(Nortmann et al. 
2018), HAT-P-11b: 0.84 A(Allart et al. 2018), HD 189733b: 0.99 A(Salz et al. 
2018), WASP 107b: 0. 84A (Allart et al. 2019; Kirk et al. 2020), HD 209458b: 
0.44 A (Alonso-Floriano et al. 2019), GJ 3470b: 1.2 A (Ninan et al. 2020; Palle 
et al. 2020). 
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Figure 9. TOI-3757 b (circled in green) in a mass—radius plane alongside other 
M dwarf planets (colored by the T.4;). We also include planets around FGK 
stars in the background, along with density contours for 0.3, 1, 3 g cm™~*. TOI- 
3757 b is the lowest-density planet orbiting an M dwarf with precise mass and 
radius measurements. 


future spectroscopic and photometric observations to place 
tighter limits on the presence of He 10830 A absorption in 
TOI-3757. 


5. Discussion 
5.1. TOI-3757 b: A Low-density Gas Giant 


TOI-3757 b is shown in Figure 9 with respect to other M 
dwarf giant planets (RK, > 4 Rs) with mass measurements at 3 
or higher. The data are taken from the NASA Exoplanet 
Archive (Akeson et al. 2013) and include recent M dwarf 
transiting planets discovered by TESS. While TOI-3757 b is 
Jovian in size (~1.05 R,), its mass is about one-quarter that of 
Jupiter (~0.26 M,); due ,” this TOI-3757 b has the lowest 
density (9 ~ 0.27 gcm °) hitherto of M dwarf planets with 
precise mass and radius measurements (>30). 

We examine different hypotheses to explain the low density 
(and puffy nature) of TOI-3757 b: 


1. Stellar insolation—The inflated radii of hot Jupiters,*’ 
ie., the so-called radius anomaly, is primarily explained 
by the high stellar insolation flux incident on these 
planets with typical equilibrium temperatures >1000 K 
(Thorngren & Fortney 2018). Here, a fraction of the 
absorbed energy from the host star is transported deep 
into the planetary atmosphere, which then through 
various mechanisms can cause its inflation.‘7 However 
as seen in Figure 10(a), the insolation flux incident on 
TOI-3757 b is not anomalously higher compared to 
similar M dwarf gas giants. Therefore, we consider it 
unlikely that the inflated nature of TOI-3757 b is due to 
stellar insolation. 


| While formally a hot Jupiter under the period based range of 1 < P< 10 
oe (Wright et al. 2012; Wang et al. 2015) TOI-3757 b does not share most of 
the planetary characteristics of hot Jupiters, with insolation flux <2 x 10° W 
m ~, equilibrium temperature <1000 K (Thorngren et al. 2016), and should 
therefore not be termed such. 

> See Fortney et al. (2021) for a review of inflationary mechanisms of hot 
Jupiters. 
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b) Planet density as a function of stellar metallicity 


Insolation (S~) 
a) Planet density as a function of insolation flux 


Figure 10. We show the planetary density of gas giants (R, > 4R,) around M dwarfs (solid colors) as a function of stellar insolation and stellar metallicity, 
respectively. In panel (a) the markers are color coded by the equilibrium temperature, whereas panel (b) is color coded by planetary radius. Additionally, TOI-3757 b 
is highlighted with a green circle. In panel (a), we show that TOI-3757 b does not receive an unusually large amount of incident flux from its host star, compared to 
other similar gas giants. In panel (b), we show that the host star has the lowest metallicity of all known gas giants around M dwarfs, which could be a potential 
explanation of its density. 


2. Low metallicity—As indicated in Figure 10(b), TOI-3757 TOI-3757 b, due to the two reasons mentioned above, we 


b has the lowest metallicity of all known gas giants 
around M dwarfs. It is important to note the caveat here 
that the M dwarf gas giants form a heterogeneous sample 
relying on different techniques for metallicity determina- 
tion,*” and hence this must be interpreted cautiously. 
While we pursue a detailed discussion on the role of 
stellar metallicity on the formation of gas giants around 
M dwarfs in an upcoming manuscript (S. Kanodia et al. 
in preparation), we discuss it briefly here. Assuming that 
the progenitor protoplanetary disk had the same metalli- 
city as the host star, the disk for TOI-3757 is ~0.3 dex 
poorer in metallicity compared to the median metallicity 
of the other transiting gas giants around M dwarfs 
(Figure 10(b)), albeit with the caveat that our metallicity 
estimate for the host star TOI-3757 has a 1o error of 0.2 
dex. Gan et al. (2022) discuss the metallicity trends seen 
in the transiting M dwarf Jovian sample, which have a 
median metallicity of ~0.3 dex, when the sample is 
updated to include the Jovian planets that have been 
confirmed since. This directly influences the formation of 
the planet in two ways—(i) There will be ~2x lesser 
material in the disk than those of comparable gas giants 
with higher metallicities (10*°? ~ 2); (ii) Yasui et al. 
(2009) suggest that low-metallicity protoplanetary disks 
have shorter lifetimes and disperse faster (than compar- 
able high-metallicity disks). This would further exacer- 
bate the problem of slow formation timescales for gas 
giants around M dwarfs, relative to the disk lifetime 
(Laughlin et al. 2004). 

The initiation of runaway gas accretion under the 
core-accretion theory requires the formation of a rocky 
core of mass ~10M, (Pollack et al. 1996) and an 
envelope comparable in mass (Meore ~ Meny), in a timely 
manner, before the disk dissipates. Under the hypothesis 
that the low stellar metallicity is driving the density of 


postulate that for such planets the process of runaway gas 
accretion did not initiate in a timely enough manner. This 
would involve a core, which though massive enough for 
runaway accretion, did not form quickly enough to 
accrete substantial gas before the disk dissipated. This 
would explain why TOI-3757 b could not form a planet 
closer to Jupiter in mass, and is a failed Jupiter. 


. Excess internal heat—We explore the possibility that 


TOI-3757 b’s radius is inflated by internal heating (as 
opposed to external insolation). Multiple studies (e.g., 
Bodenheimer et al. 2001; Burrows et al. 2007; Lopez & 
Fortney 2014; Millholland 2019) show that a planet’s 
radius is influenced by its internal heat—hotter interiors 
lead to “puffier” radii and thus lower densities compared 
to cooler interior planets of the same composition. One 
driver of a hotter interior is simply set by a planet’s age; 
younger planets naturally possess hotter interiors (Marley 
et al. 2007). The lack of a detectable stellar rotation 
period and SED fit indicate that TOI-3757 is an older star, 
and therefore we expect its core to have naturally cooled. 
If TOI-3757 b’s inflated radius is due to a hotter interior, 
it must be heated by some other mechanism. We 
investigate tidal heating created by orbital eccentricity 
as a method for inflating the radius of TOI-3757 b. While 
our eccentric orbit determination for TOI-3757 b is only 
at the ~2¢ level with an eccentricity of 0.14 + 0.06, we 
use the framework presented in Leconte et al. (2010) and 
calculate a tidal-to-irradiation luminosity ratio of ~0.05 
(assuming an eccentricity of 0.14, zero obliquity, and a 
reduced tidal quality factor—Q’—of 10°). Millholland 
et al. (2020) find that a tidal-to-irradiation luminosity 
ratio of >10~° is required for potential radius inflation 
due to eccentricity driven tides. With a tidal luminosity 
ratio of 5%, TOI-3757 b could be experiencing significant 
tidal heating from its eccentric orbit. From this power we 
calculate an interior temperature of 500 K for TOI-3757 


b. We also determine that any eccentricity >0.001 would 
lead to some tidal heating in TOI-3757 b and provide 


43 Refer to Passegger et al. (2022) for a detailed discussion on the caveats and 
challenges associated with estimating M dwarf metallicities. 
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Table 4 
List of Giant Planets with (R, > 4 R;) Around M Dwarfs, with Precise Masses (>3¢), and Spectroscopic Stellar Metallicity Estimates, as Shown in Figure 10 
Host Star Pl. Radius Pl. Mass Pl. Density Orbital Period Tege Metallicity Reference 
Rs Ms gem days K dex 

GJ 436 419+ 23.1798 1.7 3te tae 2.64388 3479180 0.1405 Turner et al. (2016) 
LP 714-47 47103 30.8113 1.636733 4.05204 395073} 0.4140-16 Dreizler et al. (2020) 
TOI-1728 5.04*9:18 26.82*3:13 1.155%0 547 3.492 3985130 0.25+¢:1 Kanodia et al. (2020) 
TOI-674 5.254017 23.6733 0.8994 9-133 1.97714 3514437 0.114007 Murgas et al. (2021) 
TOL-532 5.79* 6.18 60.38134, 1.7157*6304 2.32665 3957183 0.38 16-04 Kanodia et al. (2021) 
TOI-3629 8.297033 g2t6 0.8*9:! 3.93655 387030 0.5401 Cafias et al. (2022b) 
HATS-75 9.917913 156.051134 0.88478-68 2.78866 379072 0.5270-08 Jordan et al. (2021) 
Kepler-45 10.761133 160.49738-6 0.717 9372 2.45524 3820739 0.287014 Johnson et al. (2012) 
HATS-6 11.194031 101.04338 0.39740 083 3.32527 3724+ 18 0.279.59 Hartman et al. (2015) 
TOI-3714 11.327934 222.01108 0.857038 2.15485 3660130 0.1534 Cafias et al. (2022b) 
HATS-74 A 11.577034 464.03 1442 1.65370 -188 1.73186 3776*9 0.5179-03 Jordan et al. (2021) 
TOI-3757 12.0294 85.25785 0.2714 0.04 3.43875 3913738 0.0795 This work 


TOI-1899 12.8979 :43 209.7743333 0.547 0.08 29.02 3841432 0.317013 Cafias et al. (2020), Lin et al. (in prep.) 


Note. Data are taken from the NASA Exoplanet Archive (Akeson et al. 2013). We exclude au mic b, NGTS-1, and TOI-530 from this table due to lack of 
spectroscopic metallicity estimates. 


interior temperatures >200 K. Therefore, it is possible relationships such as Wolfgang et al. (2016), Chen & Kipping 
that tides are responsible for at least a part of TOI-3757 (2017), Ning et al. (2018), and Kanodia et al. (2019). Attempts 


b’s radius inflation. Quantifying the full impact of tidal have been made to include the orbital period (or insolation) as a 
heating on the radius inflation however requires a detailed third dimension (Ma & Ghosh 2021) in these models to capture 
investigation into the interior structure of TOI-3757 b and the radius anomaly due to host stellar insolation (Miller & 
a better estimate of eccentricity. Fortney 2011; Thorngren & Fortney 2018). However, as shown 
4. Planetary rings—Recent work by Piro & Vissapragada above, stellar insolation is not the driving force behind the low 
(2020) investigate whether rings could explain the density of TOI-3757 b, suggesting the need for higher (>3) 
apparent inflated radii of low-density superpuffs. For dimensional analysis including stellar metallicity or eccentri- 
rings to create the appearance of inflated radii due to city (among others). We also use the giant planet models from 
deeper than expected transits, the rings would need to be Fortney et al. (2007) to predict the rocky core mass for TOI- 
at an oblique angle to the planet’s orbital plane. If a planet 3757 b. While we do not have a precise age estimate, the 
becomes tidally locked to its star however, the ring models assuming a system | Gyr old, predict a core mass of ~6 
system must remain in the orbital plane of the planet (no Mz. The models for 10 Gyr (the oldest planetary population 
tilt). In this case, we would only observe the edge on modeled by Fortney et al. 2007) do not accommodate objects 
portion of the ring during transit and have a negligible as low-density as TOI-3757 b; however a planet as large as 


impact on the overall transit depth. Using Equation (11) TOI-3757 b would need to have a mass of ~140 Mz (which is 
from Piro & Vissapragada (2020), we calculate a about 1.6x the mass of TOI-3757 b) to have a core just ~1 Mz 


synchronous timescale of 5.3 Myr for TOI-3757 b. If in mass. We obtain the same result with the planetary models 
TOI-3757 b did possess rings, they would lie in its from Baraffe et al. (2008), which are also unable to 
orbital /rotation plane and not explain the planet’s inflated accommodate the low density of TOI-3757 b. 

radii. That said, we are viewing this system from a 

slightly inclined angle, which may tilt the rings into our 5.2. Atmospheric Characterization 


vantage point. We thus further investigate whether a ring 


system could remain stable given TOI-3757 b’s proxi- TOI-3757 b spans multiple unique regions of the parameter 


mity to its star. Ohta et al. (2009) note that the Poynting— space, in terms of its low bulk density and low host-star 
Robertson drag timescale is extremely short for close-in metallicity (Table 4). Understanding how TOI-3757 b formed 
planets. Using Equation (17) from Ohta et al. (2009), we and evolved could provide the context for explaining the low 


density of this unusual planet. Atmospheric characterization 
with the Hubble Space Telescope (HST), JWST, and future 
instruments could provide the necessary information to test our 
hypotheses regarding TOI-3757 b: is its low density a result of 
its formational history, or is TOI-3757 b’s radius inflated from 
present day tidal heating? 

Despite the relatively faint host star (J mag ~ 12), the ~3% 


calculate a drag timescale between 0.1 and 5 Myr for 
TOI-3757 b dependent on the assumed ring particle 
densities and sizes. Unless TOI-3757 b also possesses 
some outside method for stabilizing its ring (such as 
shepherding moons similar to Saturn), it cannot dynami- 
cally maintain a significant ring system. 


Due to the electron degeneracy pressure, Jovian-sized transit depth and large scale height of ~500 km for TOI-3757 b 
objects can span ~100x in mass, ranging from 0.3 M, (Saturn lend it a large transmission spectroscopy measurement (TSM; 
mass) to ~80 M, (brown dwarfs and very-low-mass stars). Kempton et al. 2018) of 190. As shown in Figure 11, it has one 
Therefore, 2D mass-—radius relationships are unable to of the highest TSMs of any gas giant (R, >6 Ra») with an 
accurately model Jovian-sized gas giants. The low density of equilibrium temperature <1000 K. Beyond its detectability, 
TOI-3757 b also makes it an outlier for mass—radius TOI-3757’s mass (and scale height) is constrained to >9c0 
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Figure 11. TSM of TOI-3757 b (circled in green) as a function of the planetary 
radius, alongside other M dwarf Saturn- and Jupiter-type planets (R, > 6 Re; 
colored by the equilibrium temperature). Similar planets around FGK hosts 
with equilibrium temperatures <1000 K are shown in the background. TOI- 
3757 b has the highest TSM among M dwarf gas giants (R, > 6 Re) with mass 
measurements, making it a lucrative target for atmospheric characterization, as 
discussed in Section 5.2. 
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Figure 12. We create a simulated JWST/NIRSpec-Prism transmission 
spectrum of TOI-3757 b (assuming a single transit) using uncertainties 
calculated with PandExo (Batalha et al. 2017) assuming an aerosol-free 
10xsolar-metallicity atmosphere model (blue) generated using Exo-Trans- 
mit (Kempton et al. 2017). We also plot the same underlying composition 
model plus an opaque aerosol layer at 1 mbars (green) and 0.1 mbars (orange) 
for comparison. We assume these aerosols to be gray absorbers (i.e., no 
wavelength-dependent absorption features). Based on PandExo simulations, 
we find that even with a high-altitude aerosol layer of 0.1 mbars, we should still 
detect water, methane, and carbon dioxide present in TOI-3757 b’s atmosphere 
with JWST at ~ 9o compared to a featureless flat spectrum. 


precision allowing us to resolutely analyze its spectrum and 
removing the planet’s mass as a confounding factor during 
atmospheric retrievals (Batalha et al. 2019). 

With one JWST/NIRSpec-Prism transit of TOI-3757 b, we 
should easily retrieve an abundance for water, methane, and 
carbon dioxide in its atmosphere. This combination of 
molecules provides a constraint on both the overall atmospheric 
metallicity and the C/O ratio (Moses et al. 2013; Heng 2017); 
measurements that inform where in the disk this planet formed 
(Oberg et al. 2011; Figure 12). We also show in Figure 12 that 
we should be able to observe molecular absorption features at 
wavelengths >2 j1m even with an aerosol layer at 0.1 mbars. 
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While degeneracies between aerosols and atmospheric compo- 
sitions occur, we demonstrate that TOI-3757 b is a promising 
target for future JWST observations even if it has high-altitude 
aerosols present. 

Aerosols (condensation clouds or photochemically created 
hazes) are ubiquitous across all exoplanet atmospheres; 
flattening or muting features in some while having little 
apparent effect in the atmosphere of others. By studying Titan 
in our own solar system, we know that UV flux (even a 
minimal amount) drives the production of photochemical tholin 
hazes from the photodissociation of methane (e.g., Hérst 2017). 
Even though TOI-3757 is an inactive M dwarf, it is likely still 
producing enough UV flux to create photochemical hazes in 
TOI-3757 b’s atmosphere. We would therefore hypothesize 
that TOI-3757 b has a hazy atmosphere. However, both 
Crossfield & Kreidberg (2017) and Dymont et al. (2021) found 
no Statistically significant relationship between the UV flux of 
the star and the water amplitude feature observed by HST/ 
WEC3. Alternatively, Dymont et al. (2021) do detect a 
tentative correlation between the planet’s density and the 
amplitude indicating that lower-density planets may also be 
more hazy. This trend is greatly influenced by the four 
superpuff spectra, all of which appear to be featureless (Libby- 
Roberts et al. 2021; Chachan et al. 2020; Alam et al. 2022). 
TOI-3757 b’s density is similar to these superpuffs though it is 
>8X more massive. 

TOI-3757 b is instead most similar to the low-density sub- 
Saturns HAT-P-18 b and HAT-P-19 b, both of which orbit 
early K dwarf stars (Hartman et al. 2011). Observed by WFC3, 
Tsiaras et al. (2018) retrieved a clear water abundance for 
HAT-P-18 b. However, when this spectrum was normalized by 
the scale height of the atmosphere in Dymont et al. (2021), the 
spectrum became nearly featureless. It is therefore unclear what 
we can expect for TOI-3757 b in regards to aerosol formation 
making it an interesting target for further aerosol studies. 

As discussed in Fortney et al. (2020), we may also be able to 
probe the interior temperature of TOI-3757 b by searching for 
signs of disequilibrium chemistry in the planet’s JWST 
spectrum. The model plotted in Figure 12 assumes chemical 
equilibrium; however, if TOI-3757 b possesses a hotter interior 
from tidal heating, methane will appear depleted compared to 
water, which remains unaffected. Moreover, with a strong 
vertical diffusion coefficient (K., >10*) in its atmosphere, 
ammonia would also be present in TOI-3757 b’s JWST 
spectrum. If neither methane depletion or ammonia is observed, 
this would also place an upper limit on the interior temperature 
of TOI-3757 and thus a limit on the tidal heating as a source for 
radius inflation. 

Given TOI-3757 b’s low density and extended atmosphere, 
the planet is potentially experiencing some atmospheric mass 
loss (even though TOI-3757 is fairly quiescent). Verifying this 
mass loss is occurring and quantifying the rate would enable us 
to build a picture of the planet’s atmosphere over time. The 
helium 10830 A line provides an observable mass-loss tracer. 
As we mentioned in Section 4.2, we rule out >7% excess of 
helium absorption around the planet (assuming a circular orbit), 
a feature that may be due to a lack of mass loss or a lack of UV 
stellar flux pumping helium up to its metastable state. We note 
that helium was recently detected around the low-density sub- 
Saturn HAT-P-18 b with an excess of 0.46% (Paragas et al. 
2021). Considering the similarities between the two planets, it 
is possible that TOI-3757 b maintains an exosphere beneath the 
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precision we achieved with a single transit on HPF. TOI-3757 b 
would therefore be an interesting target for future helium 
observations that achieve a higher precision measurement. 


6. Summary 


We present the discovery and confirmation of TOI-3757 b, a 
Jovian-sized planet, characterized using a combination of space 
based photometry from TESS, precise RVs from HPF and 
NEID, ground-based photometric observations from RBO, and 
speckle imaging from NESSI. With a planetary radius of 
12.0793 R and mass of 85.3*$3 Mj, TOI-3757 b has the 
lowest bulk density (9 = 0.27193 g cm) of all M dwarf gas 
giants. Additionally, its host star (TOI-3757) has the lowest 
stellar metallicity (~0.0+0.20), of all M dwarfs hosting 
transiting gas giants. We present different hypotheses to 
explain the low density of the planet, with the most plausible 
hypotheses being a formation scenario where its low metallicity 
is responsible for the delayed on-set of gaseous runaway 
accretion before the protoplanetary disk dissipated, and an 
evolution mechanism where the tidal heating causes the 
inflation of the planet due to its possibly slightly eccentric 
orbit (e ~0.14 + 0.06). 

We observed a transit of TOI-3757 b using HPF to place 
upper limits on helium 10830 A absorption by the planetary 
exosphere. Finally, we discuss how the low density (and large 
scale height) of the planet make it an excellent target for 
transmission spectroscopy. With just one transit of JWST, we 
should be able to retrieve the abundance of water and methane 
in its atmosphere, which would enable us to place limits on its 
C/O ratio. Measuring the levels of methane and ammonia 
would help constrain the interior temperature of TOI-3757, 
which can provide insight into potential tidal heating. 
Additionally, we discuss how TOI-3757 b can also be used 
to test the correlation between hazy atmospheres and planetary 
density. 
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